Tunable decoherence in the vicinity of a normal metal— superconducting junction 
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The decoherence rate of a quantum dot coupled to a fluctuating environment described by 
a normal-metal superconductor junction is considered. The density-density correlator at low fre- 
quencies constitutes the kernel which enters the calculation of the phase coherence time. The density 
fluctuations are connected to the finite frequency current-current correlations in the point contact 
via the continuity equation. Below and above the gap, at zero temperature the density correlator 
contains spatial oscillations at half of the Fermi wave length on the normal side. As the bias crosses 
the superconducting gap, the opening of new scattering channels enhances the decoherence rate 
dramatically, suggesting the possibility of tuning the decoherence rate in a controllable manner. 



I. INTRODUCTION 

The problem of decoherence in mesoscopic systems 
has remained a central issue for two decades. With the 
possible advent of quantum computers, it now occupies 
an even more important role, as the limitations associ- 
ated with various decoherence mechanisms provide the 
fundamental working limits of these devices. Propos- 
als for studying decoherence in Aharonov-Bohm (AB) 
type geometries have been madeo. Experiments in the 
last decade have provided information on the phase shift 
suffered by an electron propagating through a quantum 
doiu'cl. These interference experiments show that coher- 
ent propagation through a dot is indeed possible. A 
few years ago, decoherence was introduced artificially in 
these same AB devices by placing a quantum point con- 
tact in the vicinity of the dota. Alternatively, dephasing 
has been studied experimentally using a phase sensitive 
double dot detectora. Because quantum transport is a 
stochastic process, current noise or charge fluctuations 
act as a dissipative environment coupled to the dot, pro- 
viding decoherence without particle transfer to this envi- 
ronment. 

In the present work, the decoherence rate is computed 
for the situation where the discrete level (the dot) is 
coupled electrostatically to the fluctuating charges of a 
normal metal-superconductor junction with an arbitrary 
bias. In the Andreev regime, at zero temperature, it is 
expected that the calculation of the decoherence rate is 
similar to that of a normal metal point contact except 
that the charge of the carriers is replaced by the Cooper 
pair charged. This simple analogy fails both at finite tem- 
peratures and at voltage biases superior to the supercon- 
ducting gap. Specifically, in both the Andreev and the 
sub-gap regime, the decoherence rate depends crucially 
on the range of the potential which couples the dot to its 
environment: when this range is lowered below the Fermi 
wavelength, oscillatory terms in the density-density cor- 
relator contribute substantially to the decoherence rate. 
Even more stunning is the behavior of the decoherence 
rate as the voltage bias crosses the gap. The opening of 
new scattering channels accounts for an additional charge 
noise, and thus provokes a sharp increase in this rate. 



In principle, this could allow to tune the system from a 
"quantum" behavior to a "classical" behavior in a con- 
trollable manner. 

The computation of the decoherence rate in the present 
work, which typically involves the calculation of a 
density-density correlator in the limit of zero frequency, 
uses as a starting point known results for the finite fre- 
quency current noiseoa, rather than a direct computa- 
tion of the density fluctuations. Indeed, the charge fluc- 
tuations are directly connected to the current fluctua- 
tions using the continuity equation, which holds in sec- 
ond quantized notation, as will be shown below. 




FIG. I. A quantum dot (QD), which contains a sharp 
level, is connected to two semi-infinite leads (lower half of 
the figure) . In its vicinity, a fluctuating current flows through 
a nanoscopic (single channel), normal metal-superconductor 
junction (upper half). Broadening of the level occurs in the 
presence of an electrostatic coupling between the dot and the 
normal side of the NS wire. 

Density-density correlations have been recently com- 
puted on a formal levelcl in the geometry of Fig. ^. In con- 
trast, here one is interested both in ultra-small junctions 
and in the transition to the above-gap regime. More- 
over, it will be shown that decoherence also occurs for an 
"ideal" NS junction, the superconducting-normal metal 
analog of an adiabatic point contact. In all of the above, 
spatial oscillations of the density-density correlations are 
shown to occur, and to enter the computation of the de- 
coherence rate. To our knowledge these oscillations have 
not been described in previous work. 
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The paper is organized as follows. The model with 
its basic assumptions is described in section [□]. The 
relationship between the phase coherence time and the 
density-densitjj-iluctuations has been established by sev- 
eral authorscJ~li3 and is reviewed briefly in Appendix |A|. 
The main emphasis is put on the analysis of the density- 
density correlations of the normal-metal superconduct- 
ing quantum dot, hrst in the Andreev regime (section 
IV). Expressions which apply to the case where the dis- 
order potential is smooth - a superconducting adiabatic 
point contact - are provided. Finally, the case of a bias 
which is superior to the gap is exposed in section |y|, and 
its consequence on the decoherence rate is illustrated nu- 
merically using the Blonder Tinkham Klapwijk (BTK) 
modelEl 



II. DEPHASING MECHANISM 
AND NS JUNCTION 

A (small) quantum dot, which for practical purposes 
here is represented by a discrete, sharp energy level, is 
located in the proximity of a (single channel) normal 
metal-superconductor (NS) junction (Fig. Q). The dot is 
connected to two semi-infinite leads, and can in principle 
be pwA of an interferometer such as the ones studied in 
RefsJ-ra. The dot is coupled by Coulomb forces to the 
fluctuating charges located in a NS point contact (SPC). 
To be specific, it is assumed that the electrostatic cou- 
pling is restricted to the normal side of the junction, as 
the dephasing mechanism is expected to be more efficient 
when the bias voltage on the NS wire is below the super- 
conducting gap. 

The system is described by a Hamiltonian which char- 
acterized the dot, the NS junction and the coupling be- 
tween the two: 



e'e 



dxU{x)^ (x)ip(x) 



(1) 



where i^ix) an( i c (V' ( x ) an d are fermion creation and 
annihilation operators in the NS junction and in the dot. 
The potential U(x) depends on the location of the charges 
in the NS junction. While it originates from long range 
Coulomb forces, in practical situations it is screened by 
the surrounding metallic gates. Later on, different ranges 
will be specified for U (x) in order to observe their conse- 
quence on the decoherence r ate . . 

The "standard" proceduraI2rlI2l for computing the de- 
coherence rate is to identify an exponential decay in 
time of the dot electron Greens function. Even though 
there is no electron leak from the dot to the NS junction, 
the level in the dot acquires a finite width due to the 
coupling with the fluctuations in the junction. The dot 
Greens function thus acquires a non-oscillatory compo- 
nent: 



where is the decoherence time, and the bracket no- 
tation ( )tvs implies that an average over the NS envi- 
ronment has been taken. The following result for the 
decoherence rate neglects the back-effect of the dot on 
the NS junction: 
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- = ^ I dtK(t) 



(3) 



with the Kernel defined as: 



K(t) = -J dxx J dx 2 U(x 1 )U(x 2 ) 

x ((p(x 1 ,t)p(x 2 ,0) + p(x 2 ,0)p{x 1 ,t))) (4) 

where the notation (( )) means that the average densities 
have been subtracted. 

The main issue to compute this rate is to specify the 
density-density correlator using the property that elec- 
trons and holes on the normal side, and possibly quasi- 
particles on the superconducting side, are scattered at 
the NS junction. The scattering matrix specifies the am- 
plitudes of the transmitted/reflected particles. It allows 
to give the asymptotic behavior of the elecitpn and hole 
wave functions away from .the junctionErllJIla. Here, the 
notations of previous worktj are used for convenience. 

On the normal side of the NS point contact, the elec- 
tron and hole wave functions associated with a particle 
a = e, h which originates from side j = N, S are given 
by: 



5 ah 

(hv + )^ SNjeh( 



Vj, a (x, t) 



ryyj(e + s Njhh e ) + 



(to_)!/2 



{hv_y/i SN]hee 



ik—x 



(•5) 



(6) 



where the electron and hole momenta are specified by 
k± = ^/2m (ps i e)/h, with p,$ the chemical potential 
of the superconductor, which assumed to be large com- 
pared to both the gap A and the bias eV. Sij a p is the 
amplitude for a particle of type [3 = e,h which is in- 
cident from side j to be scattered as a particle of type 
a = e, h in reservoir i. When discussing the Andreev 
regime, both electrons and holes are incident only from 
the normal side, and the indices i, j — N will therefore 



G(t) = (T[c(t)J(0)}) NS oc e-*/ r *e- ieot ^ 



(2) 



be dropped in sections [II and IV, but will be restored in 
section [y| 

Finally, the electron/hole distribution function on the 
normal side are given by: /jv e ,ft = {1 + ex p[( e T 
eV)/kB®}}~ 1 while on the superconducting side, inci- 
dent holes and electrons have the same distribution func- 
tion, with V = 0. 
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III. CURRENT AND DENSITY FLUCTUATIONS 

As described before^, the statistical average of the 
current-current operator is-iobtained by performing the 
Bogolubov transformation!!!!] on the current operator. In 
the past, it was for the most part computed at equal lo- 
cations in connection with the current noise across the 
NS junction. For the decoherence time, it is necessary 
to keep the two locations xi and x*. separate. However, 
as the current correlators in Refs.LrO are expressed with 
separate time arguments t\ and ti , the generalization to 
separate spatial arguments is straightforward: 

2*2 /. + 00 .+<» 

= ^Y,J o I d£d£ '{ 

Q, 

*[(«8 du* & ) tl , xi + {vp dv%)t U xt] 
x[(u*p du & )t 2 , X2 + («| dVa)t 2 ,x 2 ] 

x{ u *pdv* & ) tl , xi [{ u &dvp)t2,x2 + {u^dVa)t2,x 2 ] 

H±-umi-fp(^y {E+E ' Kt2 - ti)/h 

(7) 

where a = (a, i) is a short hand notation combining 
the reservoir from which the particle, electron or hole, 
is incident. Eq. (^) constitutes the starting point for 
computing both finite frequency noise and the zero fre- 
quency noise in the presence of a local harmonic pertur- 
bation, such as in the Non Stationary AB effect in NS 
junctionsQ which was recently detected experimentallyo. 
Here it constitutes the starting point for the computation 
of the density-density correlator, and it is valid both in 
the Andreev regime and above gap, provided that the 
proper distribution functions are specified on the super- 
conducting side. 

The continuity equation allows to relate the current 
operator to the density operator as it holds in second 
quantized form: 

p(x, oj) = — V.J(x, w) (8) 

This allows to write a connection formula between the 
nonlocal current noise correlator and the density-density 
correlator at finite frequency: 

f +oc dt 

({p(x 1 ,u>)p(x 2 ,-u>))) = —e^ 

J-oo ^ 

x d Xl d X2 ((I(Q)I(t))) (9) 



Here the u> — density-density correlations will be 
needed. Taking the derivative with respect to the po- 
sitions, the to 2 term in the denominator is canceled, for 
all bias regimes, giving a finite contribution to the density 
fluctuations. 



IV. DENSITY CORRELATOR IN 
THE ANDREEV REGIME 

In the Andreev regime, there are only two types of par- 
ticles (electrons and holes), so the current-current corre- 
lator at finite temperatures and bias takes the form: 

e 2 h 2 f +oc f +OQ r 
{(I(x 1 ,t 1 )I(x 2 ,t 2 )))A = iz-2 / dede'\ 

Lm Jo JO 

/ e (e)(l-/,(e'))e 4(£ '- £)(t2 - tl)/?1 
X[(u h du* e ) tuXl + (v h dv*) tuXl ] 
x [« du e ) t2 , X2 + {v* h dv e )t 2 ,x 2 ] 

X[(u e dU* h ) tl ,xi + («e dV* h ) H , Xl ] 

x[(u* e du h ) t2 , X2 + (v* dv h ) t2 , X2 ] 
+/ e (e)/, i (e')e- l(£+e ' )(t2 - tl)/?1 

x(u* h dv*)t 1:Xl [(u e dvh)t 2 ,x 2 + {uhdv e )t 2 ,x 2 } 
+/, l (e)/ e (e')e- l(£+£ ' )(t2 - tl)/?l 

x(u*dvl)t uXl [(u h dv e ) t2 , X2 + {u e dvh)t 2 ,x 2 \ 
+(l-/ e (e))(l-A(e'))e l(e+£ ' )(t2 ^ l)/?l 

x[{u e dvh)t x ,xx + (u h dv e )t 1 ,x 1 }(u* h dv*)t 2 ,x 2 
+(l-/, l (e))(l-/ e ( e '))e 4(£+£ ' )(t2 -* l)/?l 

X[{u h dv e )t u x 1 + (UedVh^.xAiKdvDt^^j 

(10) 

To proceed, one makes use of the unitarity of the scat- 
tering matrix, combined with the time reversal symme- 
try of electrons and holes in the Andreev regime (energy 
dependence of scattering coefficients neglected): s^ e = 
— s e h and s* e = Shh- Taking the Fourier transform of 
the current-current correlator, the following relation be- 
tween the incoming and the scattered wave numbers is 
obtained: 




Moreover, the standard simplifications are made (ps S> 
lo): once the spatial derivatives are taken and the u) 2 pro- 
portionality in d xl d X2 {{I{xi,ti)I(x2,t2))) is identified, 
the wave vectors of electrons and holes are assumed to 
be equal to fci?. 
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The density correlator in the Andreev regime becomes: 

9 p+oo , 



((p{x 1 ,uj)p(x2,-ui)))a = 



f Jo 



f e (e)(l-f h (e-u))Q(e-cj) 
x4|s eft | 2 [|s ee | 2 + s* e e 



xA\s eh \ 2 [\s ee \ 2 + s* ee e 2lk ^ +s ee e 
-f e {e)f h {-e + u)Q(-e + u) 

X\s eh \ 2 {\ - (\s eh \ 2 - S ee Shh) + See 

+f h (e)f e (-e + u)G(-e + u) 



^ — 2ik F x 1 _|_ ^2ik F (x 2 — xi)j 
-2ik F x 2 _|_ g2ik F (xi— x 2 )l 



— 2ik F X2 



Shhe 



2| See | 2 {|s ee | 2 + 1 + s ee (e- 2ik ^ + i/ 2e - 2ik ^) 
+s* ee {e 2lkFX1 + l/2e 2tkFX2 ) + cos(2k F (x 2 - a*))} 



+s* e e 2lkFX2 {l + s* ee e 2lkFX1 ) 



-2ikp r -Ci 



(i 



(12) 



Further assuming that the temperature fcg6 < A, and 
neglecting the energy dependence of the scattering ma- 
trix coefficients, the thermal integrations are performed 
in Appendix [b|. 

At zero temperature (uj > 0) the only interval which 
survives is the one specified by the Fermi functions 
/ e (e)(l - f h (e - lj)) (first term in Eq. (|l|)). At low 
frequencies, and fcs0 = 0, the double integral in energy 
is eV and the correlator becomes: 



the different role played by charge and current fluctua- 
tions for the decoherence rate. This constitutes the N-S 
analog of the adiabatic point contact studied in Ref.tZl. 
This situation assumes a sharp NS interface with per- 
fect Andreev reflection, adjacent to a scalar potential on 
the normal side which varies slowly on the scale of the 
Fermi wave length . While the Andreev reflection pro- 
cesses have a unit probability, it is necessary to take into 
account the dependence of the electron and hole wave 
function amplitude on the normal side due to the smooth 
disorder potential following a quasi-classical/ WKB ap- 
proximation. For simplicity, here only results at zero 
temperature, in the Andreev regime are presented. 

Instead of redefining the electron and hole wave func- 
tion in this limit, it is more convenient to directly sub- 
stitute the expressions for the semi-classical matrix ele- 
ments in the density correlator of Eq. (12): |sft, e | 2 = 1 
and See = Shh = 0. The wave numbers of electrons and 
holes, which now have acquired a spatial dependence, still 
satisfy the relationships of Eq. ( |ll| ) . 

For the SAPC, the only contribution which survives 
is the one which is directly proportional to the Andreev 
reflection probability: 



((p{x 1 ,u))p(x 2 ,-uj))) 



I i KB 



16e 2 f d J e {e){\- f h {e-Lo)) 
h 2 



y/v+(x 1 )V-(x 2 ) 



exp 



HI + k(x))dx 



({p(xi,u)p(x2,-uj)))a,u=o = 



2e 3 V\s eh \< 



VW + v+{xi)){v'_(x 2 ) + v-{x 2 )) 



(15) 



X IM 2 + s* ee e 2ikpX2 + s ee e~ 2lkFXl + e 24A *<*2-*i)] 

(13) 

Note the remarkable fact that the low frequency 
density-density correlator has an oscillatory spatial de- 
pendence with wavelength X F /2. These oscillatory terms 
give a significant contribution to the kernel of Eq. (jij) 
when the envelope function U(x) is short ranged. To- 
gether with the above approximations, the decoherence 
rate in the Andreev regime can be expressed in terms of 
the Fourier components U(q) of the envelope potential at 
q = and at q — ±2k F : 



Specifying that the distribution functions are step 
functions, this simplifies to: 



{{p(x 1 ,u)p{x 2 ,-u))) 



16e 3 F 



WKB 



exp 



h 2 yjv+{xi)v-{x 2 ) 

Q i(k'(x) + k{x))dx 



VW^i) + y+(xi))(v'_(x 2 ) + v_(x 2 )) 



(16) 



— r = lim 



tfvj? Jo 



tfc/ e (e)(l-/fc(e-w))e(e-w) 



x[s* he (s' hh U(0)+^U(-2k F )) 

-<h(4eU(0) + U(-2k F ))] 
x[ She (s h * h U(0) + U(2k F )) 

- s '; h ( See u(o) + u(2k F ))}, 



(14) 



where prime denotes quantities evaluated at e — to. 
The case of a superconducting adiabatic point con- 
tact (SAPC) is now considered momentarily, to highlight 



Note the analogy of Eq. (|lj) with the expressions de- 
rived for normal, adiabatic point contaculj. This illus- 
trates that even for a NS junction with ideal transmis- 
sion, the decoherence rate does not vanish, as substan- 
tial density fluctuations are present, although the current 
fluctuations are reduced due to the Pauli principle (which 
operates on electrons and holes on the normal side). 



V. DECOHERENCE RATE AT 
ARBITRARY BIAS 

In the previous expressions, contributions where quasi- 
particles are transmitted into the superconductor were 
discarded. Here, the calculation of the density-density 
correlator proceeds as before, choosing for simplicity zero 
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temperature, a constraint which excludes some combina- 
tions of fNe.h and fse,h- The calculation proceeds in a 
similar way as in the Andreev limit, except that processes 
involving quasi-particle emission from the superconduc- 
tor now contribute. Operating the simplifications on the 
wave vectors as before, and using the continuity equation, 
the density correlator becomes: 

((p(x 1 ,oj)p(x 2 , -u))) - ({p(x 1 ,u)p(x 2 , -u))) A = 

2e 2 ,+oo 



de\ 



h 2 V 2 F j 

f Ne (e)(l- fse(e-Lo))Q(e-uj) 

X ( s NShe s *NNhe + s NSee( e 2lkpXl + S NNee )) 
X ( s NShe s NNhe + S NSeei^"^ + S NNee )) 

+f N e(e)(l- fsh(e-u))e(e-uj) 

x ( s NShh s *NNhe ~ s NSeh{ e 2lkpXl + S* NNee )) 
X { s NShh s NNhe - S^ Seh (e 2lkFX2 + SNNee)) 

+.fNe{e)fse{-£ + w)9(-e + UJ) 

XS NSee S *NNhe( S NSee S NNhe + SNShe(e 2lkpX2 + s' NNee )) 

+fNe{e)fsh(-e + w)9(-e + w) 

XS NSeh s *NShe( s NSeh s NNhe + S N Shh{e 2lkFX2 + s'iViVee))} 

(17) 

where ((p(xi,uj)p(x2,—^)))a is the density-density cor- 
relator with Andreev scattering contributions only (yet 
the energy integral ranges from to eV, above the gap). 

Taking into account that the decoherence rate in- 
volves only zero frequency density-density correlations, 
and that quasi-particle energies are always positive, the 
decoherence rate reads: 



f 



1 



t£ n 3 v 2 F jo 

(sNSheS NN heU(0) 

+ S NS ee{U{-2k F ) + S* NNee U{0))) 
x {s*NShe s NNheU(Q) 

+ s NSee 

(U{2k F ) + S N NeeU{Q))) 
+ (sNShhS N Nhe U ( Q ) 

-s NSeh {U{-2k F ) + s* NNee U(0))) 

x ( s *NShh s NNheU(0) - 

s NSeh (U(2k F )+s NNee U(0)))} 



de\ 



(18) 



In order to enquire about the effect of the oscillatory 
terms, a Gaussian potential profile with a specific width 
£ is chosen: 



U(x) = j-e & 
U(k) = -%e~^ 



2ir 



(19) 



where £ represents the screening length of the Coulomb 
interaction due to the surrounding metallic gates. For 
nanoscopic dots and junctions, or alternatively for the 
large wavelengths which apply to semiconductors-2D 
electron gas structures, it is becoming conceivable that 
X F could become larger than £. 

With this particular choice, the decoherence rate can 
be expressed as: 







= * I 
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ri) 


2irh 3 v 2 F Jo 



-8tt 



2_e^ 

{sNSheS* NNhe + S N See(e ^ + S* NNee )) 

X ( s *NShe s NNhe + S* N See( e * F + SNNee)) 

-Stt 2 



+ (sNShhS NNhe — SNSeh(e 



- 2 A 



S NNee)) 



X ( s *NShh s NNhe ~ S* NSeh {e A -f + S N Ne 



»} 

(20) 



where the Andreev decoherence rate reads: 



t£U 2 2nh 3 v 2 F jo 



eV 



del 



[s*NNhe( s NNhh + e X f ) - S NNeh (s* NNhe + 6 V )] 

X [sNNhe(s*NNhh + e F ) ~ S* NNeh (s N Nee + e A F )] j 

(21) 



The results are now illustrated by plotting the de- 
coherence rate as a function of bias voltage, for differ- 
ent values of the potential range £. Because first, it 
is relevant to enquire about the role of disorder in the 
NS junction, and second, the complete energy depen- 
dence of the scattering matrix coefficients is required, 
a model with a minimal set of parameters, the BTK 
modeiLj, is chosen. Expressions for the scattering ma- 
trix elements of this model are knowrfi there, the same 
model was chosen to enquire about the singularities in 
the frequency dependent noise in a NS junction, for ar- 
bitrary biases. Recall that the NS junction is described 
by a delta function barrier V(x) = VbS(x) and a step- 
wise pair potential A(x) — A0(x), both located at the 
normal metal-superconductor interface. The strength 
of the (normal) barrier is represented by the variable 
Z = mVBjh 2 k Fl which for intermediate values between 
high and low transparencies is of order 1. 
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FIG. 2. Decoherence rate £ 2 [/?>L?o] _1 as a function of bias 
voltage in the Andreev and sub-gap regime, for barrier trans- 
parencies: Z = 0.1 (dashed line), at £ = 0.15Af (top) and at 
\ = IOAf (bottom); Z = 1.0 (full line) at £ = 10A F (top) and 
at £ = 0.15Af (bottom). 

For high to intermediate barrier transparencies, the de- 
coherence rate is plotted in Fig. |[ This rate is normal- 
ized to the interaction strength squared ([/o/£) 2 , which 
allows to plot the curves corresponding to Z = 0.1 and 
Z = 1.0 on the same scale. In addition, results are il- 
lustrated for two values of the range £ for comparison: 
£ = IOAf corresponds to a range where the negative ex- 
ponentials in Eqs. ( ppf ) and ( |2l| ) can be neglected; on 
the opposite, for £ = 0.15Af these exponentials give a 
significant contribution. 

Deep within the Andreev regime eV <C A, the deco- 
herence rate have a linear dependence on the bias (not 
shown), as the scattering matrix coefficients are essen- 
tially energy independent. The rates computed with 
these parameters show no noticeable change in slope once 
the voltage bias crosses the gap. This linear dependence 
is expected for large, sub-gap biases, because transport 
across the NS junction then becomes dominated by sin- 
gle quasi-particle transfer: similar behavior was recently 
observed for the finite frequency noise of NS junctional. 

It is difficult to predict the behavior of -rj 1 when the 
energy dependence of the scattering matrix coefficients 
becomes noticeable, i.e. for biases near the supercon- 
ducting gap. In Fig. ||, no qualitative changes of the 
different decoherence rates are observed until eV ~ A, 
above which the rates show a gradual crossover to the 
linear dependence on the bias. Yet, the results plotted 
in Fig. H show that the relative magnitude of these (nor- 
malized) decoherence rates depends in a non trivial man- 
ner on the potential range £: for intermediate barriers 
(Z = 1.0) the oscillatory terms in the density-density 
correlator (Eq. |j]) tend to reduce £ 2 [ r 0^o] -1 , while for 
high transparency barriers, they increase the decoherence 
rate. In addition, for eV ~ 1.9 A the curves for £ = IOAf 
corresponding to intermediate and weak barriers cross, 
and the rate corresponding to Z = 1.0 is larger beyond 
this. 



A 1.5A 

FIG. 3. Same as FigJ^ for opaque barriers (Z = 10.): 
£ = IOAf (dashed line), and £ = 0.15Af (full line). 

Most interesting is the behavior of the decoherence rate 
for opaque barriers (Z = 10.), which is illustrated in Fig. 
||. For both potential ranges £ = IOAf and £ = 0.15Ap, 
t^ 1 displays a clear cusp as the bias crosses the gap. This 
cusp corresponds to the opening of new scattering chan- 
nels - electron and hole quasi-particle transfer from the 
normal metal to the superconductor - which gave previ- 
ously a gradual change - and no obvious crossover region 
- for higher transparencies (Fig. ||): here, the crossover 
region is confined in a very small interval with a width 
lesser than a few percent of A. Above the gap, the linear 
dependence on the bias is recovered. 

Note that the decoherence rate should in principle be 
highest when the current undergoes strong temporal fluc- 
tuations: one naively expectsEj that it is related to the 
shot-noise fluctuations^ ~ T{\ — T). Here, T stands 
either for the Andreev reflection probability (below the 
gap) or for the transmission probability of quasi-particles 
(sub-gap regime). Indeed, the curves corresponding to 
Z = 10 (Fig. H) could not be plotted on the same scale 
as in Fig. @ because the corresponding rate is too small. 

The drastic change depicted in Fig. || could possibly 
be exploited to switch the decoherence rate from "large" 
to "small" , thus allowing to control the degree of coher- 
ence - from quantum to "classical" (phase incoherent) 
transport - in the neighboring quantum dot (Fig. |l|). 
In addition, note that for existing NS junctions opaque 
barriers are likely to be the norm rather than the excep- 
tion, which renders an experimental check of this effect 
plausible. 



VI. CONCLUSION 

The low frequency density-density correlations of a NS 
point contact have been extracted from the finite fre- 
quency noise characteristic, using the continuity equa- 
tion. This density correlator enters the computation of 



G 



the dephasing rate of an isolated quantum system with 
a discrete spectrum. 

In the shot noise regime, the density fluctuation func- 
tion contains terms which oscillate with half a Fermi 
wavelength. Because this length scale is considerably 
larger in semiconductors than in metals, and given the 
ongoing progress in nano-fabrication techniques, it is 
suggested that the detection of such oscillations, or their 
effects in the dephasing rate, could indeed be possible. 
The role of these oscillatory terms is noticeable for weak 
disorder even in the (sub-gap) Andreev regime. 

Expressions of the density-density correlator described 
in section IV for a NS junction were generalized to the 



case of a Superconducting Adiabatic Point Contact, us- 
ing a semi-classical scheme. Note that this implies that 
decoherence can even occur when the NS junction located 
near the quantum dot has ideal charge transmission prop- 
erties (pure Andreev scattering). 

While the computation of the decoherence rate beats 
similarities with that of a normal-normal point contactlia, 
no simple extension so far has been provided when 
both Cooper pairs and quasi-particles are transmitted in 
the superconductor. The competition between the two 
charge transfer processes has been a central issue in the 
present work. 

Regardless of the disorder strength, the decoherence 
rate rises sharply as the voltage bias is increased above 
the superconducting gap. Yet for opaque barriers, this 
transition is even more dramatic, as indicated by the cusp 
obtained in Fig. [| This crossover is due to the opening 
of additional scattering channels, which renders the fluc- 
tuating environment more noisy. It is conceivable that 
this increase of the decoherence rate could have some 
potential applications. 

In particular consider a quantum bit (qubit), imple- 
mented by a quantum dot apparatus: the coupling to a 
neighboring NS junction can render the system classical 
(besides the tunneling processes in and out of the dot) 
or quantum in a perfectly controllable manner. Decoher- 
ence, reduced abruptly by lowering the voltage bias be- 
low the gap, could be used to determine when the quan- 
tum evolution of the qubit is supposed to start, or to 
end. Other applications include the classical resetting of 
a quantum computer. 

The present approach has dealt with a single channel 
NS wire. Clearly a generalization to a many channel wire 
would bring this proposal closer to experimental realiza- 
tions, in carbon nanotubes and so on. Yet, the main goal 
here has been to demonstrate the variability of the de- 
coherence rate on the bias voltage as it crosses the gap. 
Because a specific model (BTK) was used to plot the de- 
coherence rate, the additional complexity of dealing with 
more than one channel could have rendered these results 
less explicit, but extensions are indeed possible. 

Finally, for more immediate experimental applications, 
it is necessary to work with superconducting materials 
which have a small enough gap that imposing a voltage 
bias at the NS junction does not cause substantial heat- 



ing on the normal side. Indeed, the present theoretical 
approach assumes that the transport across the junction 
is purely elastic, in the mesoscopic regime. Such super- 
conductors with gaps of the order of hundreds of mK 
are readily available, and allow presently to observe a 
crossover behavior in the current noise above and below 
the gapEI 
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APPENDIX A: CALCULATION OF THE 
DECOHERENCE RATE 

The starting point for the decoherence rate is the dot 
Green function: 



G(t) = -i<r[c(t) c t(o)]> 



(Al) 



Operators are in the Heisenberg picture. Specifying the 
time evolution: 



< c{t)J(0) >= 

< e -^ot/h Tte --k J* dt' j dxU(x)^ (x,t)i/i(x,t) C /Q\ c t /Q\ 



> 
(A2) 



Following Levinsonli!li3, it is assumed that while the 
NS point contact has an effect on the dot, the reverse is 
not taken into account. The Greens function acquires a 
non-oscillatory time dependence: 



< c(t)c f (0) >= e- leat exp[-(j){t)} 



(A3) 



where the decoupling between the dot and point contact 
degrees of freedom allows to compute the average as in a 
Gaussian process: 



4>{t) ~ - In 



<T t e 



% J* dt' J dxU{x)ip^ (x,t)tp(x,t) 



>NS 



\ j dt 1 j dt"K{t' - t") 
Ti Jo Jo 



(A4) 



The kernel K is computed assuming that the quantum 
dot does not perturb the NS point contact: 

1 



K(t) = -J dx x J dx 2 U(x 1 )U(x 2 ) 

x ((p(xi,t)p(x 2 ,0) + p(x 2 , 0)p(X!,t))) 

(A5) 

This expression connects the decoherence rate Eq. (|^) to 
the density fluctuations in the superconducting SPC. 

The long time be hav ior of the non oscillatory time de- 
pendence in Eq. (A3) is obtained by considering the 
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density-density kernel. The kernel is characterized by a 
time scale t c which identifies for which times the corre- 
lations still survive. Therefore, in the limit t — > oo the 
second integral over t" will saturate to a constant : 



i 

2ft 2 
t 

T 4> ' 



d£ 



dt"K{t") 



(A6) 



assuming that for long times, one can replace the second 
integral over the whole time domain. 



2|s ee | 2 {M 2 + 1 + s ee (e- 2lfcFXl + l/2e" 

- 

) 



+s* ee (e 2lkFXl + l/2e 2lkFX2 ) + cos{2k F (x 2 - a*))} 



-2ikpx-2 



(i 



-2ikpx 



■)]} 



(B3) 



At zero temperature this result corresponds to Eq. 
(ED- 



APPENDIX B: FINITE TEMPERATURES 



In this appendix, the energy integrals which enter in 
Eq. ( |T^ ) are computed analytically. 

In the Andreev regime , the elements of the scattering 
matrix are assumed to be weakly dependent on e, so one 
only needs to compute the integrals of the Fermi Dirac 
distributions, which have the general form: 



de- 



1 



1 + e P(cTeV) 1 + e -f3(e-u±eV) 

1 f l + e^ eV ) 



(Bl) 



de 



1 



1 



1 + e P(^eV) I + e l3(-e+uj±eV) 



—M- 



1 + e±? eV 



(3{e^ - 1) 1 1 + eP(" ±eV ) J 

(B2) 

Combining the thermal factors, the density-density 
correlator can be computed as : 



((p(ll,u)p(l 2 ,-w))) 



h 2 h 2 v 2 F 

1 + ef eV 



l /3(e^(«- 2«V) _ l) m \ 1 + e /3(^-eV) J 

x[|s ee | 2 + s* ee e 2ikpX2 + s ee e- 2lfcFa;i + e «*H*a-*i)] 
1 + -i 



-rit»(- 

1 ll 



2 j_ c * p 2ik F xi j_ ^ e ~2ik F x 2 _j_ e 2ife_p (rci -a: 2 ) 



X[|S ee | +S ee e 



1 f 1 + eP eV 1 

+ fl( e ^ - 1) I 1 + e /3("+eV) J 



x{l - (\s eh \ 2 - s eeShh ) + s ee e~ 2 * k ^ + s hh e 2lk ^} 



— — 

1 U 



/3|seh| 2 (e^-l) ll + e/ 3 ^- el/ ) 
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